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ABSTRACT 

Motivated by the anomalous dust-to-metal ratios (DTTVI) derived in the literature for -y-ray burst (GRB) damped Lya absorbers 
(DLAs), we measure these ratios using the dust-depletion pattern observed in UV/optical afterglow spectra, associated with the 
interstellar medium (ISM) at the GRB host-galaxy redshifts. Our sample consists of 20 GRB absorbers and a comparison sample 
of 72 QSO-DLAs, overall at redshift 1.2 < z < 4.0 and down to Z = 0.002 Zq metallicities. The DTM in QSO- and GRB-DLAs 
increases both with metallicity and metal column density, spanning ~10-1 10% the Galactic value and pointing to a non-universal dust- 
to-metal ratio. The low values of TfTM. indicate that low metallicity systems have lower dust fractions than typical spiral galaxies 
and perhaps that the dust in these systems is produced inefficiently, i.e. by formation in asymptotic giant branch (AGB) winds and 
by grain growth in the low-metallicity regime, with little contribution from supemovae (SNe). On the other hand, some GRB- and 
QSO-DLAs show high TfTM values out to z ~ 4, requiring rapid dust-production, such as in SN ejecta, but also in AGB winds and 
via grain growth for the highest metallicity systems. GRB-DLAs overall follow the TfTM. properties of QSO-DLAs, GRBs probing 
up to higher column and volume densities. For comparison, the TfTM. that we derive for the SMC and LMC are ~82-100% and ~98% 
of the Galactic value, respectively. The literature TfTM. value of the low-metallicity galaxy I Zw 18 (< 37%) is consistent with the 
TfTM distribution that we find. The dust extinction Ay increases steeply with the column density of iron in dust, A'(Fe)j„st, calculated 
from relative metal abundances, confirming that dust extinction is mostly occurring in the host galaxy ISM. Most GRB-DLAs display 
log A'(Fe)d„s, > 14.7, above which several QSO-DLAs reveal molecular hydrogen, making GRB-DLAs promising candidates for 
molecular research. 

Key words. Gamma-ray burst: general - ISM: abundances - (ISM:) dust, extinction - (Galaxies:) quasars: absorption lines 

1. Introduction Large Magellanic Clouds (LMC) have Ay/A?(//)x[SMClI] ~ 

^ ■ 2.8 X 10-" mag cm^ and Ay/A?(//)x[LMC]i ~ 5.3 x lO^'^ mag 

. : The mterplay between the production/recycling of metals and ^m^, corresponding to ~ 50% and ~ 95% of the Galactic 

in the formation/evolution of cosmic dust is strongly connected y^j^g^ respectively. Such measuremen ts are tvpicallv verv dif- 

^ with star formation (e.g., |Gall et al.||201l4 [ Cortese et al. | |2012|) ^^^^^ for distant galaxies. iDai & Kochane^ (l2009l) 

^ , and has a fundamental role for the chemical enrichment of the studied the dust-to-metal ratio in gravitationally lensed galaxies 

. . ■ interstellar medium (ISM). at z < 1 by comparing the extinction curve and the X-ray ab- 

> ' One way of estimating the dust content al ong the li ne of sorption in the spectra of multiply-imaged background quasars 

sight is from the optical dust extinction Ay (e.g.,|Pel|l99l. The (QSOs). They find the dust-to-metal ratio to range from ~ 10% 

low-ionization gas content is typically traced by the neutral hy- to a few times the Galactic value, but with large uncertainties. 

; drogen A^(Hi), while a possibl e probe of the metals is the soft Dust-to-metal ratios increasing with metalh city were recently 

- - X-ray absorption A^(//)x (e.g. [Morrison & McCammo^|1983|: derived for super-solar metalUcitv galaxies bv lBrinchmann et al.l 

[Wilms et alj|200g). Ay/NiH i) can then be used to infer the dust- ([2oTl . base d on observed CO absorption and the theoretical ex- 

to-gas ratio. However, this ratio does not provide a natural way pectations of lCharlot & Longhettil (|2001l) . 

to quantify the amount of dust in a galaxy's interstellar medium Long-duration (> 2 s) GRBs can offer the unique possibiUty 

as It does not include the metals in the gas. On the other hand, directly probing the gas, dust and metal content in their star- 

the dust-to-metal ratio, i.e. the fraction of metals in dust, repre- foming host galaxy, out to very high redshifts. Their afterglows 

sents the ability of a system to produce or destroy dust grains, ^^^^^ bright enough to allow absorption-line spectroscopy 

Besides, it can be derived regardless of the H i constraint, which jhat c an reveal the H i gas and metal column densiti es in the ISM 

IS often not available (see Sect[2|. We refer to DTM to indicate jSavagho et alJ [2001 IVreeswiik et aP [lool iBergeret al.1 
the dust-to-metal ratios normalized by the Galactic value. 

The dust-to-metal ratio in the Galaxy is Av/N(H)x[Gal] ■ Derived from the observed E(B-V)/N(m from lMartin et al.l (fl989h 
5.6 X 10 mag cm^ ([Predehl & Schmitt 1995), in agree- and iDraind (12003!) and convert ing it to N(H)x/Av assuming Ry = 

3li 



ment with the me asurements based on UV absorption fines 2.9. 3.2 i lPeilU992l : lDrainell2003l) and a mean [Zn/H] = -0.64,-0.55 
dBofilin etalJll978l) . Tfie SmaU Magefianic Cloud (SMC) and for SMC and LMC respectively jWeltv e"taDll997l . ll999h . 
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2006t iProchaska et"an 120071: iLedoux et al] 120091: iFvnbo et al.l 
20091). In add ition, the well-understood power law continuum 



(e.g.. ISari et al . 1998; Piran 1999) facilitates the determination 
of the intrins ic Ay and N(H)x- The latter is measured for most 
Swift b ursts ([Evans et al"] 2p09[) and Ay can b e studied out to 
z > 6 (ISchadv et al.l 120101: IZafar et al.ll201 Ibh . In contrast, the 
same measurements are more difficult to make towards QSOs, 
due to the presence of broad emission lines contaminating a large 
fraction of their optical/ultraviolet continuum. The samples of 
GRB absorbers that are typically followed up for spectroscopic 
studies are biased against the most dusty systems that have been 
dimmed by large Ay's or faint GRBs. A fraction o f GRB host 
galax i es is massive and dusty ( Kriihler et alJ 1201 it iHunt et aT] 
120111: iRossi et all 120121: iPerlev et all l2013l) . Nevertheless, the 
current (biased) samples is still representative of the majority 
of the GRB host population, i.e. l ow-mass, low-metalU city and 
actively star-forming galaxies (e.g. lSavaglio et al ] l2009l) . 

IZafar et al.1 (1201 lal) and ISchadv et al"] (1201 Oh investigated 
the dust- to-metal ratios for a large sample of distant GRB- 
selected galaxies using Av/N{H)x- They found much lower 
Av/N{H)x in GRBs than what had been observed in the Local 
Group, up to two orders o f magnitudes, confirming the results of 
iGalama & WiiersI (l2001h . This implies that either GRB-DLAs 
have a different dust-to-metal ratios or that Ay and N{H)x probe 
different regions along the line of sight. While the dust is ex- 
pected to lie in the cold low-ionization medium, the X-ray ab- 
sorption traces the metals regardless of their ionization state. The 
evidence of ionization induced by the GRB on the surrounding 
medium (in the immediate environment of the GRB for most 
bursts and out to hundreds of pc for absorbing sy stems with low 
Cia et al ] l2012t IVreeswiik et al.ll20T3h . confirms that 
indeed the line of sight analysis can be affected by ionization and 
dust destruction. Thus, the Av/N{H)x may not provide the best 
estimate of the dust-to-metal ratio for GRB host galaxies. Recent 
estimates based on line of sight Ay, A^(Hi) and the metallicity 
suggest that the dust-to-metal ratios in GRB an d QSO absorbers 
are ro ughly constant around the Galactic value (IZafar & WatsonI 
I20TI . 

An alternative approach to calculate the extragalactic dust- 
to-metal ratios is based on the dust-depletion pattern revealed 
by UV/optical absorption-line spectroscopy. Using such meth- 
ods, the dust-to-metal ratios in QSO-DLAs and a few GRB- 
DLAs were foun d to be consistent with the Galactic value 
or slightly lower (Vladilo"1998'; 'Savaglid l200Tl ISavaglio et al] 
12003 : Savaglio & Fall 2004; Vladilo 200#. One of the advan- 
tages of using optical absorption lines rather than dust extinction 
to determine the dust content is that the absorption lines refer to 
a particular redshift and thus they are not integrated along the 
whole line of sight. This is particularly crucial for low metallic- 
ity systems with an intrinsically low dust content and for high-z 
targets, whose lines of sight likely cross intervening galaxies, 
each contributing to Ay- 

In this paper we develop a new dust-depletion-based method 
to derive the dust-to-metal ratios, cross-checked with the 
ISavagliol (l200ll SOI hereafter) method, and apply it to a sam- 
ple of 20 GRB absorbers (observed at medium to high spectral 
resolution) and 72 QSO-DLAs (all with high-resolution spec- 
troscopy). This approach solely relies on the observation of op- 
tical/UV low-ionization absorption lines associated with DLAs 
and therefore there is no or little ambiguity on the location or 
ionization state of the absorbing gas/dust: we purely focus on 
the bulk of the ISM in the host galaxies, which is mostly far- 
away from the GRB (e.g., hundreds of parsecs from the burst 



IVreeswiik et al.ll2007l 1201 ih and therefore not affected by ion- 
ization due to the GRB or the progenitor. 

The paper is organized as follows. Our sample is defined in 
Sect. |2] and the methods of our analysis are explained in Sect. 
|3] The results are presented in Sect. |4]and discussed in Sect. |5] 
Finally, we summarize our conclusions in Sect. |6] Throughout 
the paper we adopt ions cm"^ as the linear unit of column den- 
sities N. We refer to relative abundances of two chemical ele- 
ments X and Y defined as [X/Y] s log ^ - log ^|f^. For the 
reference solar abundances appearing in the second term of this 
formula, we use the valu es of Asplun d et al. (2009) and follow 
the recommendations of iLodders et al.l (l2009l) by adopting the 
photospheric estimates for the more volatile elements, the mete- 
oritic estimates for the less refractory elements , or the average 
between them (see details in lDe Cia et alj|20T2l) . 



2. The GRB-DLA sample 

We select all GRB absorbers with a constrained estimate of the 
column density of Fe 11 and of at least one ion among Zn 11 and 
Si II from the largest to date sampl e of G RB absorber column 
densities collected bv ISchadv et al.l (1201 ll) . The resulting sam- 
ple is composed of 20 GRB absorbers in the redshift range 
1.2 < z < 4.0 and spanning a wide range of metallicities 
(from Z ~ 2Zq down to Z ~ 0.005 Zq), listed in Table [T] 
Most absorbers have an H i estimate from Lya, all b ut one with 
log A^ (H i) > 20.3 and therefore classified as DLAs (Wolfe et al.l 
mm . One exception is GRB 090323 with logA^(Hi) = 19.6 for 
one of the two absorbing systems that composes its line pro- 
file. The column density estimates of the GRB sample are de- 
rived from absorption-line spectroscopy, using curve-of- growth 
an alysis for the spect ra with low or mid resolution (see details 
in Is chad V et al.ll20Tll and references therein) a nd Voigt-profile 
fitting for a handful of hi gh-resolution spectra (IVreeswiik et al.l 
l2007tlLedoux et allllool . 



3. Revisiting tlie dust-to-metal ratio for GRB-DLAs: 
the iSavaqlioit2001|) recipe and a new method. 

The dust-to-metal ratio can be estimated from the depletion in 
dust of different elements, observed with optical/UV absorption- 
line spectroscopy, in particular from Fe 11 and Znii lines. While 
H I is fundamental to derive the overall metallicity, it is not nec- 
essary to derive this ratio. Fe 11 and Zn 11 are two elements that 
have a similar nucleosynthetic history but very different deple- 
tion properties. Fe is typically heavily depleted onto dust grains , 
while Zn is only mil dly depleted (e.g.. lSavage & Sembachll996l : 
iLedoux et al.l |2002|) . Zn thereby being a metallicity indicator. 
Thus, these are the most suitable elements for this analysifl Si 
and S can be useful dust-depletion indicators as well, bu t should 
be tre ated with caution, especially in the case of S (e.g.. lJenkinsl 
l2009l) . as discussed below. 

3.1. The S01 method 

The D7M (referred to as k/kEin SOI) can be derived from the 
depletion pattern, i.e. from the observed abundance of different 
elements, as compared to the depletion pattern observed in dif- 
ferent Galactic environments: halo (H), dis kn-halo (DH), warm 
disk (WD) and cool disk (CD), as listed in lSavage & SembachI 



^ Cr could be used as well, but it is not as easily observed as Fe. 
^ did, in lSavaglio et alj Hmk and lSavaglio & Falll (|2004^ 
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Table 1. dust-to-metal ratios normalized by the Galactic value (DTTA) derived with the method described in this paper and 
SavagUo's (2001, SOI). 
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0.18 
0.21 
0.17 
0.07 
0.04 
0.03 
0.27 
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0.05 
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0.27 
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D 
H 
D 
DW 
DW 
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H 

HWC 
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DW 
all 
D 
C 
all 
all 
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D 
H 
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Notes. Uncertainties refer to a 1 cr confidence level, while 3 cr limits are reported for unconstrained measurements. " Observed abundance of the 
element X relative to solar. * Abundances corrected for dust depletion. When available, [Zn/H](ot is considered to be the metallicity of the gas. 
Column density of iron in dust-phase (see text). Metallicity as fitted with the SOI method (Eq.[T]l, when H but not Zn is co nstrained. " Galactic 
enviro nments fitted with the SOI method: halo (H), disk+halo (D), warm disk (W) and cool disk (C), as characterized bv ISavage & Sembachl 
( fT99^ . When more than one environment is possible, 2)7yVl(S01) and logZ/Zg are averaged among the allowed values. ^ Dilferent absorbing 
systems towards the same GRB are reported separately. 

References. All column densities, z and Ay are taken from the compilation bv lSchadv et al.l j201 Ih . 



(Il996h . Following the recommendation of iJenkinsI (l2009l) . we 
sligh tly update the mean depletion values listed in table 6 of 
ISavag e & Sembac^ ( Il996h by adding -0.03, -0.04 and -0.01 
to the values for Fe, Si and S respectively, here and through- 
out the paper The abundances with respect to H or Zn observed 
in GRB-DLAs are fitted simultaneously assuming the Galactic 
depletion patterns, providing a best-fit £)77V( and environment 
type. When both H i and Zn 11 are constrained we treat the metal- 
licity [Zn/H] as an input. Instead, when H i is available but not 
Zn, we fit for the metallicity Z/Zq as well. In practice, we iter- 
ate among the different Galactic environments / and sum over 
different metals X to minimize the total difference between the 
observed abundances and the reference abundances 

[X/H]rf = log(l + DTM ■ (lO'^'j - 1)) + log(Z/Zo), (1) 

where dX j is the depletion of a metal X f or a given environmen t 
J. More details can be found in SOI and lSavaglio et al.l (l2003h . 
We exclude upper and lower limits on the abundances from the 
fit. We include measurements for at least two metals X among 
Zn, Si and Fe. The results of the depletion pattern fit for each ab- 
sorber are shown in the Appendix (Fig. IA.ll and rA.2l for the case 
with and without H i constraints respectively). Sulphur is not in- 
cluded in the fit because it has shown controversial depletion 
patterns in different Galactic environments (Jenkins 2009) and it 
shows an erratic depletion behaviour in QSO-DLAs as well. The 
inclusion of silicon must also be taken with caution, as the intrin- 



sic abundance of Si with respect to Fe and Zn may not entirely 
depend on dust, but instead it could in principle be boosted by a- 
element enhancements in star-forming galaxies. This leads to a 
possible degeneracy between the best-fit DTM, metallicity and 
Galactic environments. We explore this possibility of an intrin- 
sically lower Si by calculating the DTMs for slightly lower Si 
column densities than observed (log A^(Si) - 0.2 dex). The effect 
of such change can be mostly accounted by ;) setting an upper 
error on the DTTVt of at least 0.2 for best-fit WD and DH envi- 
ronments, when Si is included in the fit and /;) setting a lower 
error on the DTM of at least 0.2 when Zn is not constrained. 



3.2. A new method to derive the dust-to-metal ratio for DLAs 



Here we derive the absolute dust-to-metal ratio (dtm) based on 
fundamental definitions and applying depletion corrections to 
the observed relative abundance of Fe with respect to Zn. The 
basic idea starts from the definition of the observed [Fe/Zn], 
where the observed gas-phase column density of an element X is 
N{X)obs = A^(-'^)tot - A^(-'^)dust- If Zn is assumed not to be depleted 
into dust and if Fe and Zn share the same nucleosythesis (total 
of gas- and dust-phase [Fe/Zn]tot - 0), then a simplified version 
of the dust-to-metal ratio follows c/fmsimpiified = 1 - 10'^^/^"l"''^ 
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In order to include the possible depletion of ZiJ3, we assume 
that the depletion of an element in dust is proportional to a deple- 
tion strength factor F*. This factor represents the overall strength 
of the dust depletion and it is unique to each line of sig ht, similar 
to what has been found for the Galaxy (|jenkinsll2009l e.g., their 
fig. 7). In the Galaxy, systems with little dust have g - and 
almost no depletion of Zn, but the dustier systems F«,g - 1 and 
even Zn is significantly depleted. We note that the Ft scale for 
DLAs does not necessarily have to be the same as in the Galaxy, 
because DLAs typically have average lower metallicity and pre- 
sumably lower dust contents. These linear depletion sequences 
reflect the tendency of each element to condensate in dust grains, 
i.e. the y are related to the conden sation temperature of each el- 
ement ( lFieldlll974l: lJenkinsll2009l) . Given that the slopes (c/Fe, 
dZn, dSi et cetera) of these linear relations are specific to each 
element, the relative composition of the dust grains changes de- 
pending on how dusty is the environment. 

While the slopes and offsets of these linea r depletion se - 
quences have been well measured for the Galaxy (|Jenkinsl2009h . 
we don't know in principle how these relations change for lower 
metallicity environments such as most DLAs. The observed off- 
sets in the Galaxy are interpreted as pre-existing dust grains that 
cause some le vel of depletion even for the less dusty systems 
(|jenkinsll2009l) . For the GRB and QSO sample of absorbers, we 
assume no offsets on the depletion sequences {6x = at F« = 0), 
meaning that for DLAs we assume no pre-existing grains. This 
assumption is justified given the relatively high redshift and the 
low metallicity of QSO- and GRB-DLAs, only a handful be- 
ing above [Zn/H] ~ -0.3 in our DLA sample. Furthermore, 
the evidence that most metal-p oor DLAs sh ow [Zn/Fe] ~ 
(iNoterdaeme et al.i JQQ8; Mo laro et al.l l2000l) strengthens this 
hypothesis. 

In this way, the depletion of Zn and Fe can be described as 

6zn = [Zn/H]obs - [Zn/H] intrinsic = dZn ■ (2a) 
^Fe = [Fe/H]obs - [Fe/H] intrinsic = dFe ■ F^, (2b) 

where [Zn/H] intrinsic and [Fe/H] intrinsic are the total abundances in 
the gas and dust phases, and when F« = 1, 6zn - dZn. Since 
the observed (depleted) abundances are associated with the gas- 
phase metals, the fraction of an element X in gas fx is 



fx = 10-* 



Figure [T] visualizes these relations and the quantities we derive 
below. Note that we don't assume solar abundances or fixed 
grain compositions, but rather rely on the observed relative abun- 
dances [Fe/Zn]obs as a function of the depletion strength. Since 

F, = (^Fe - Szn)KdFe - dZn) = [Fe/ZnUJ(dFc - dZn), (4) 

the dust to metals ratio dtm is then defined as 

dtm = !^^E£>±1 = 1 _ lO'*^^ = 1 - 10[F'=/Z"]ob.^Fe/(^Fe-^Zn) (5) 

A^(Fe)tot 

This definition of dtm assumes only two linear dust depletion se- 
quences for Fe and Zn respectively, and that Zn and Fe have the 
same nucleosynthetic history ([Fe/Zn]tot = 0), so any deviation 
of [Fe/Zn] from solar is due to dust depletion. 

The slopes of the depletion sequences correspond to the ex- 
pected level of depletion in the most dusty {Ft = 1) DLAs sys- 
tems. We assume dFe = -0.95, dSi - -0.26 and liZn = -0.11, 



Some level of Zn depletion should be taken into account, given that 
GRBs are likely to reside in star-foitnin g regions, some of them with a 
potentially significant dust content (e.g.. lPerlev et al.l[2013h . 




Fig. 1. A scheme displaying two linear depletion sequences of 
Zn and Fe with the depletion strength factor Ft. Despite Ft is not 
known for a given line of sight, it can be traced by the observed 
[Zn/Fe] ratio. 

i.e. th e average depletion along d iskn-ha lo lines of sighfl of F e 
and Si lSavage & Sembachl(ll996h and Zn lRofli & Bladeslll995h . 
This translates to ~ 90% of Fe, ~ 45% of Si and ~ 20% of Zn be- 
ing depleted into dust for the most dusty DLAs. The assumption 
of steeper slopes, s uch as the Gala ctic ^Fe = -1.28, c/Si = -1.14 
and ^fZn - -0.61 (|Jenkinsll2009h . would overestimate the dust- 
corrected metallicities (see Sect. l5.1l for a discussion on the de- 
pendence of our results on the assumed slopes). Thus, we limit 
our analysis to the slope values mentioned above and include 
an uncertainty of 0. 1 on the slopes in order to allow for some 
freedom. 

Finally, we normalize the dtm values with respect to the 
Galactic value 

DTM = dtm/dtm{G), (6) 

where dtm{G) - 0.89 (see the Appendix lAl for its derivation). 

Another useful quantity that can be derived from the above 
formalism is the column density of iron in dust, A^(Fe)dust (col- 
umn 9 of TablelU, which provides an estimate of the dust content 
in the ISM, independently from Ay. Given that 



N{FsU = A'(Zn)t< 



N{Fe) 



(3) and 



N(Zn\ 



[NiZn))^ 
A?(Zn)obs N(ZnU, 



20[Zn/Fe]obs«'Zn/(dZn-dFe)' 

the dust-phase column density of iron is 
A^(Fe)dust = A^(Fe)tot • dtm = 

_ jQ[Fe/Zn]„fo _ ^-|^Q[Fe/Zn]„(,j rfFe/Ci/Zn-f/Fe) _ j-j 

'A^(Fe)^ 



(7) 



(8) 



■ Ar(Zn)o 



A^(Zn) 



(9) 



4. Results 

4. 1 . Comparison among DTM methods 

The DTM derived with the new and the SOI methods for our 
GRB-DLA sample are listed in Table[T]and plotted against each 
other in Fig.|2] for comparison. The two methods show a general 
good agreement and thus we consider them as equivalent (as also 



^ Disk-Hhalo lines of sight: HD 18100, HD 167756. 
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Fig. 2. A comparison between the UTM derived with the new 
and the SOI methods. The symbols indicate which metal was 
used as a reference to calculate the [X/Fe] ratio, where [Zn/Fe] 
is the most reliable ratio because it is less sensitive to intrinsic 
abundance changes. The filled symbols highlight the DLAs with 
high-resolution spectroscopy, which typically have more robust 
constraints. The vicinity to the equality (dashed) line shows the 
agreement between the two estimates. 



confirmed by a Spearman rank tes10). The results of this test are 
listed in Table |2l There is a possible small deviation towards 
higher £)77V((S01) above 1, for three GRB absorbers, namely 
GBR 010222, GRB 050401 and GRB 090323, absorbing system 
2. We discuss whether such deviation could affect our results in 
Sect. |5] The overall reliability of both methods is strengthened 
by the agreement of the two different TfTM derivations, because 
they are independent from each other 

Beside the UTM ratios, the SOI method also provides esti- 
mates of the metallicity, in case the H column density is con- 
strained but not Zn, and the environment ty pe (halo, warm disk, 
cool disk, disk + halo, as characterized by iSavage & Sembach 
Il996h . These are listed in the last two columns of Table [T] Note 
that when Si is constrained but not Zn, the metallicity fitted from 
the depletion pattern is always higher than the observed [Si/H], 
as expected from the depletion of Si. Among the best-fit environ- 
ments, a cool disk is somewhat disfavoured for GRB absorbers, 
although a single best-fit environment could not be constrained 
for about half of our sample. 

In the early stages of the current analysis we also derive d 
the dust-to-metal ratio using the prescription of IVladilol(ll998h . 
but we do not include those results because of the assumption of 
a constant grain composition (fixed fraction of each element in 
dust). Nevertheless it is worth reporting that those dust-to-metal 
ratios are consistent with the results presented here. 



4.2. Comparison with QSO-DLAs 

In order to compare the GRB- and QSO-DLA populations 
we calculate the DTM for a sample of 72 QSO-DLAs from 
Noterdaeme et al. (2008), all of which have been observed with 
the high-resolution Very Large Telesco pe (VLT) Ultraviole t and 
Visual Echelle Spectrograph (UVES, iDekkeret al.ll2000l) . We 
choose to compute the £)77Vl of QSO-DLAs with our new 
method, because it is easily derived from the [X/Fe] reported in 
Noterdaeme et al. (2008). The overall sample of GRB and QSO- 
DLAs spans a wide range of metallicities, from Z ~ 2 down 
to Z ~ 0.002 Zq. The average observed [Zn/Fe] in the QSO and 
GRB samples are 0.4 and 0.9, respectively, with standard devia- 
tion 0.2 and 0.5. This suggests that GRB environments are over- 
all more dusty than QSO-DLAs. Indeed, only less th an ~ 10% of 
QSO- DLAs show significant reddening due to dust dKhare et al.l 
I2012h . 

Figure [3] shows the distribution of DTyVIs, derived with the 
new method, in GRB-DLAs (black) and QSO-DLAs (grey) as 
a function of metallicity (left panel) and the Zn column density 
(right panel), corrected for dust depletion using Eq. |8] The deple- 
tion corrections are typically very small, up to ~0.3 dex. In Fig. 
|4]we also include those targets that have measurements of the 
Si column density but not of the Zn column. The DTM that has 
been estimated from Zn should be considered more reliable, be- 
cause the [Zn/Fe] is less aff'ected by intrinsic changes of relative 
abundances than [Si/Fe]. 



4.3. Comparison witli fhe Magellanic Clouds 

We calculate the £)7yV( values for two lines of sight through the 
SMC and the LMC, using the Zn and Fe column densities de- 
rived from the absorption lines observed in the literature. For the 
SMC, we co nsider to the absor ption in spectra of the background 
star S k 108 dWeltv et al.lll"997lf l Sk 155 and Sk 78 dWeltv et al.l 
1200 Ih. For the LM C we use the absoiption towards SN 1987 A 
(IWeltv et al.llI999P l While the choice of a few specific line of 
sight per system may not represent the overall properties of the 
whole cloud, we note that typically only one line of-sight can 
be observed towards both GRB and QSO as well. The derived 
DTM are -82-100% and -98% of the Galactic value for the 
SMC and the LMC, respectively, as displayed in Figs. [3] and |4] 
Given the similar [Zn/H] and logA?(Zn) of the SMC and LMC 
absorbers, it is interesting to note that the LMC has a higher 
[Zn/Fe], DTM, and ultimately a larger amount of dust with re- 
spect to SMC, although these diff'erences are not 3 cr significant. 



5. Discussion 

5.1. Strengths, caveats and implications of our DTM method 

The estimates of the dust-to-metal ratio presented in this paper 
rely purely on the metal absorption lines observed in the op- 
tical/UV afterglow spectra. These are selected to be at the GRB 
redshift, typically arising in th e ISM of the host galaxy, hun dreds 
of parsecs from the burst (e.g. IVreeswiik et al.ll2007[ |20I ll) . 

The main assumption of our D77A method is the existence 
of linear depletion sequences in DLAs, similar to what has been 



* This indicates only how tight is a monotonic correlation between 
two quantities, but they do not need to be necessarily linear. 



' We refer to the mean observed [Zn/Fe] = 0.50 ±12 and [Zn/H] = 
-0.64^551^ as sociated wit h the SMC, consistent with [Zn/Fe] = 0.66±20 
measured bv lSofia et al.l (120061) for the same line of sight. 

** We refer to the mean [Zn/Fe] = O.W° \l and [Zn/H]= -0.55 ±0.08, 
averaged among the velocity components associated with the LMC. 
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observed in the Galaxy (|Jenkinsll2009l) . These sequences de- 
pend on the condensation temperature of each chemical element 
and naturally explain different dust-grain compositions in dif- 
ferent environments. While t he slopes of th e Galactic depletion 
sequences are well studied (iJenkinsI |2009|) . the corresponding 
slopes in DLAs are unknown. Here we discuss how much our 
results depend on the assumption of these slopes. 

At the low metallicity end, i.e. in the least dusty DLAs, 
the low depletions are not very sensitive to changes in slopes 
(see Fig. [TJ, while the zero offsets of the linear sequences are 
fully justified (see I3.2l i. For the most dusty DLA systems, the 
dependence on the slopes of the depletion sequences can be 
more significant. However, by comparing with the Galactic val- 
ues, we find that the possible range of slopes is limited. The 
slopes of the depletion sequences for DLAs must be shallower 
than in the Galaxy to provide reasonable dust-corrected metal- 
licities in the DLAs. Indeed, using the Galactic values of the 
depletion-sequences slopes would give a general small increase 
of TfThX (up to ~ 0.15), but imply excessive metallicities. In 
the most extreme cases, the Galactic slopes would lead to dust- 
corrected metallicities of up to Z ~ VZ© and Z ~ I2OZ0 
for GRB 000926 and GRB 090926A (which are individual ab- 
sorbing systems with [Zn/Fe]- and [Si/Fe]- based TfThX respec- 
tively), which is unrealistic. These numbers are even more ex- 
treme for more dusty absorbers such as towards GRB 090323, 
where the Galactic slopes lead to dust-corrected Z ~ 30 Z© and 
Z ~ 10"^ Zq for the two separate absorbing systems. These results 
are not justifiable, regardless of the degeneracy of the H i con- 
tent between the two absorbing systems. On the other hand, the 
slopes that we assumed for DLAs provide dust-corrected metal- 
licities Z ~ 1.1 Z0, Z ~ 0.01 Z0, Z ~ 3.4 Zq and Z ~ 2.5 Zq for 
the systems mentioned above, respectively. Thus, the depletion 
of Si, Fe and Zn are likely milder in DLAs than in the Galaxy. 
The £)77V1 values given in Table [T] are probably reliable within 
the errors, despite the uncertainties on the slopes of the depletion 
sequences in DLAs. 

Interestingly, the depletion of Fe and Zn of the most dusty 
DLAs (with depletion factor ~ \ and above) are similar to 
what has been observed in the least dusty lines of sight in th e 
Galaxy (with Galactic depletion factor F»,g ~ 0, Jenkins '2009^, 
such a s disk-Hhalo lines of sight studied by fSavage & SembachI 
(Il996h . One intriguing possibility is that the depletion factor 
could evolve continuously from the least depleted dust-free 
DLAs to the most dusty Galactic lines of sight, where the de- 
pletion sequences would start at = with no offset and mild 
slopes and gradually evolve into the observed Galactic depletion 
sequences. 

5.2. The mM distribution 

The TfTM distribution of the GRB-DLAs in our sample in- 
creases with metallicity and spans ~10%-1 10% of the Galactic 
value. In addition, the observed trend points to low dust frac- 
tions for low-metallicity systems, not consistent with a constant 
Galactic £)TM at an y metallicity as recently suggested (e.g., 
\Zafar£rWatw^\20lA) . These results confirms that Av/N(H)x 
estimates lead to inconsistent dust-to-metal ratios in GRB-DLAs. 

The comparison between GRB- and QSO-DLAs in Fig. [3] 
reveals a continuous distribution of the DTM between the two 
populations, indicating an overall similarity of the ISM. The 
mean D7M for GRBs and QSOs are 0.8 and 0.6 respectively, 
with standard deviation 0.3 and 0.2. Given that GRB lines of 
sight preferentially cross through inner (and more dusty) regions 
of their star-forming galaxies, while QSO lines of sight typically 



intercept the haloes of in tervening galaxies (e.g. lProchaska et aP 
120071: iFvnbo et alJl2009l) . GRBs can probe the DTM distribu- 
tion up to one order of magnitude higher column densities than 
QSO lines of sight. In the Galaxy, the average volume den- 
sity of the gas is tightly correlate d with the dep letion factor 
(F^ c = 0.772 + 0461 ■ log(«(H)) [jeiikinsl 12001 . For DLAs, 
we can derive the average volume density from the observed 
the depletion factor Ft (Eq. |4|, using the relation observed in 
the Galaxy. While the actual values of («(H)) that we find for 
DLAs may be unreliable, because the relation might evolve with 
metallicity and be different for DLAs, the comparison between 
the densities in QSO- and GRB-DLAs can be meaningful. While 
QSO-DLAs have 0.001 < <n(H)) < 10 cm-\ several GRB- 
DLAs exceed 10 cm"^, confirming that GRBs are typically lo- 
cated in denser regions than QSO absorbers. 

The DTM increases with metallicity (Fig.|3]and|4] left pan- 
els). Despite the large scatter, the increasing trend is clear up to 
[X/H] ~ -I. Above this metallicity the distribution may be con- 
sistent with a roughly c onstant dust-to-metal ra tio for DLAs, as 
found - based on - b\ jZafar & W atsonV2013!). At lower metal- 
licities we find low DTM values in both QSO- and GRB-DLAs, 
clearly not constant and not consistent with the Local Group, 
pointing to a non-universal dust-to-metal ratio. 

Th e discrepancy between our results and the Ay-derived 
mM (IZafar& Watsonll2013[ , but with one single constrained 
datapoint at very low metallicities) might arise from an over- 
estimate of the intrinsic dust content when estimated from the 
line of sight Ay of the lower metallicity and less dusty systems. 
For these low-dust systems the additional extinction due to fore- 
ground absorbers along the whole line of sight could start to 
contribute. Indeed, the observed foreground E(B - V) towards 
QSOs typically ranges between ~0.001 and ~0.1 |Menard et al] 
I2OO8I) . i.e. 0.003 < Ay < 0.3, so this additional source of dust 
extinction can be significant when compared to the intrinsically 
low Ay of very low metallicity systems. On the other hand, our 
dust-depletion analysis does not suffer from contribution of in- 
tervening systems and therefore provides reliable estimates of 
the dust content of low metallicity systems. 

In Figs. [3] and |4] we include a literature estimate of the 
mM of the blue compact dwarf galaxy I Zw 18 - one of the 
lowest-metallicity systems known i n the local universe ([O/H] = 
-1.52, Mdust/Me < 8.1 x 10'^: iHerrera-Camus et all 1201 2|) . 
where we consider the dust- to-metal mass ratio Mdust/^met = 
10i°g('Wd„„/MH)-[O/H] jjjjjj normaUze it by the Galactic value 
(Mdust/MH)Gai = 0.0073 (iDraine et al.ll2007h . The mM < 0.37 
of I Zw 18 is consistent with the trend that we find for GRB- 
and QSO-DLAs, supporting the idea that low metallicity sys- 
tems have lower dust fractions than typical spiral galaxies. 

The mM distribution increases with the metal column den- 
sity as well (represented by the Zn or the equivalent Zn column 
density, right panel of Fig. 13) and flattens out above log A^(Zn) ~ 
13. One issue that needs to be clarified is whether this flatten- 
ing is physical or have been artificially introduced by our anal- 
ysis. As noticed in Sect. |4l three mM values might be slightly 
underestimated by our method (see Fig. |2]i, all of which have 
logA^(Zn) > 13. However, given the magnitude of the effect 
these three points are unlikely to prove crucial. The observed 
flattening of the 2)7A1 at higher metal columns might suggest 
that the availability of metals in the ISM can drive dust forma- 
tion/growth up to some stable level where metals are largely 
accessible, as for the Galaxy. Dust-destruction due to GRBs 
might also play a role in limiting the amo unt of dust (e.g., 
IWaxman & Drainell2000l: fDraine & Hadl2002l) . 
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[Zn/H]joT log A^(ZnII)joT 

Fig. 3. Left: The distribution of the UTM with dust-co rrected TZn/HlxoT (the in trinsic metalUcity), for the GRB-DLAs (black) in 
our sample and the U VES QSO-DLAs (gray) sample of iNoterdaeme et al.l ( i2008h . The MW, SMC (3 lines of sight) and LMC TfTM 
values are displayed for comparison, where the latter two are two representative values derived from the observed [Zn/Fe] along 
two lines of sight (see text for details). Filled symbols highlight the measurements dra wn from high-resolution spe ctroscopy. The 
thick green arrow shows the literature value (derived from the Mdust/A^H measured bv iHerrera-Camus et al.ll20T2l) for the nearby 
blue compact dwarf galaxy I Zw 18. Right: The distribution of the TfTM with the total column density of Znn, corrected for dust 
depletion. 




-2-10 1 11 12 13 14 

[Z/H]tot log7V(ZnII)TOT 

Fig. 4. Same as Fig. [3] but including the TfTM values that have been estimated from Si and thus are less reliable. When Zn 
is not directly constrained we display the equivalent Zn ii column density derived from the Si ii column, assuming solar relative 
abundances. The symbols indicate which metals are used in [X/Fe] to derive the UTM. 
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Table 2. Spearman's rank correlation coefficients p and signifi- 
cance. 
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Fig. 5. The TfTM. evolution with redshift. The symbols are the 
same as in Fig.|4l 



The significance of the correlations described above is con- 
firmed by a Spearman rank test, as listed in Table |2] indicat- 
ing that DTTVt is most tightly correlated with [Zn/H]. There is a 
much weaker correlation between the TfTM distribution and the 
Hi column density, as shown in Fig. lA.3l in the Appendix. This 
confirms that the dust traces the metals rather than the total-gas 
content. 

While the UTM values in QSOs are uniformly distributed, 
we cannot exclude the presence of two populations of GRB hosts 
divided by their UTM properties, i.e. a standard one following 
the trends that have been observed for QSO-DLAs and a sec- 
ond group of GRB-DLAs with a lower UTM. This is most no- 
table in the D77V1 distribution against the metal column den- 
sity, although the gap is less evident when including TfTM esti- 
mates derived from Si as well. It is difficult to assess whether 
this separation is a real effect, given the small-number statis- 
tics. If true, these low-D7yVt GRB-DLAs (mostly GRB 071031 
and GRB 0804 13 A) may correspond to systems having a smaller 
amount of dust for their metal content, possibly due to sputter- 
ing or sublimation of dust grains in actively sta r-forming regions 
(iTielens et alJll994tlWaxman & Drainell2000l) . Further observa- 
tions are clearly needed to investigate this issue. 

We do not observe any trend of the IXTM with the mass, star- 
formation rate and specific star-formation rate of the host galax- 
ies, although the sample for which all these estimat e s are avail 



able is limited to a handfu l (IVreeswiik et al.ll2004t IChen et al 
2009 1: ISavagho et alJ 12009 1 : iKriihler et alJ l201 it ISavaglio et al 
2012). 



5.3. On the origin of dust in DLAs 

We further explore the evolution of the DTM with redshift, dis- 
played in Fig.|5] Combined with the metallicity information, this 
provides some further clues on the dust production in DLAs. 



5.3.1. At low redshift 

All the estimated UTM values for GRB-DLAs below z ~ 1.7 
(6 out of 20) are high, around the Galactic value. While the 
metallicity of these systems has not been measured from ab- 
sorption lines because Ly-a is not redshifted into the observ- 
able window, the Zn ii column densities of these absorbers are all 
above log A^(Zn ii) > 1 3 (with the excep tion of GRB 080330 with 
logA^(Znii) ~ 12.8 iD'EliaeTall 120091) . Low-z GRB absorbers 
follow the UTM distribution with the metal column discussed 
above and it is reasonable to expect that they follow the UTM 
distribution with metallicity as well. If so, this would indicate 
that at redshift below z < 1 .7 most GRBs occur in galaxies with 
metallicities Z > O.IZq, possibly because the mean metallicity 
of star-forming galaxies has increased with cosmic evolution. 
However, some caution should be exerted because the number 
of observations at low redshift is limited. 



5.3.2. At high redshift 

While above z > 1 .7 the UTM distribution shows a large scatter, 
we confirm the evidence for some high-UTM g alaxies to lie at 
high redshift, such as GRB 090323 at z = 3.6 dSavaglio et al 
2012) or QSO-DLA Q 1441+2737 at z = 4.2 t^edoux et aT 
200®. The presence of dusty systems at these redshifts suggests 
that a rapid dust-production mechanism may be required for the 
dust to have formed well within a Gyr One such fast mecha- 
nism is the condensation of dust grains in cooled SN ejecta (e.g., 
(Morgan & Edmunds 2003t), where the contribution of SNe on 
dust production i s dominant at earl y times (for time-scales be- 
low ~ 200 Myrs, iGall et al.|[201 lah . AGB stars start to recycle 
their metals into the ISM produc ing silicates grains after ~ 40 
Myrs (iDi Criscienzo et a"iJl2013h and carbon-ric h dust after a 
few h undreds Myrs of main-sequence evolution (iGalliano et al.l 
I2OO8I) . Grain growth in the ISM is typically considered a slow 
process to produce dust, but its time-scale decreases with in- 
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creasing metallicity (the growth time-s cale t is of the order o f 
~60 Myrs for solar-metalUcity systems. [Hirashita & Kuoll201 ih . 
Thus, dust production via grain growth cannot be excluded for 
those systems at high redshift with relatively high metallicity, 
such as GRB 090323. While the current observation of the UTM 
at z ~ 4 cannot put stringent limits on the origin of this dust, 
future observations of D77V1 at redshifts above 5 will be cru- 
cial to assess the extent o f SN contribution at high redshift 
dMorgan & Edmuridsll2003l) . The method described in this pa- 
per in principle can allow the UThX to be measured up to very 
high redshift. The observation of rapidly-produced dust (within 
30-170 Myrs) in z > 6 QSOs suggest that SNe dominate the 
dust producti on at these redshif ts, with a negligible contribution 
of AGB stars (iGall et al.ll20TTbl) . 

5.3.3. At low metallicities 

The low TfTM. of Z < O.IZ0 DLAs impHes that these systems 
have a lower fraction of dust with respect to the higher metallic- 
ity systems. This could be due to a lower efficiency of dust pro- 
duction and/or a higher efficiency of dust destruction. Below Z ~ 
O.IZ0 grain growth starts to be less efficient and its time-scale 
becomes larger (r of hundreds of Myrs to G yr, depending on 
the m etallicity and the grain-size distribution; iHirashita & Kuol 
I2OI ih . The efficiency of dust-production from AGB stars is 
also lower in low-metallicity systems, b ecoming a negligi- 
ble s ource of dust below Z ~ O.OOSZ© (iDi Criscienzo et all 
I2013|k while SNe dust production is less sensitive to metallic- 
itv dGaUianoet al..,2008l) . 

One possible cause of a higher efficiency in dust destruc- 
tion could be a harder and more intense radiation field pro- 
duced by enhanced massive-star formation in low metallicity 
environments, such as for GRB host galaxies for instance (e.g. 
ISavagUo et ani2009l) . However. lOall et al.l (1201 1 ah showed that 
simulated dust masses are higher for initial mass functions with 
more high-mass stars, meaning that SNe are overall more effi- 
cient at producing than destroying dust. An alternative cause of 
a higher level of dust destruction could be a grain-size distribu- 
tion skewed towards sm aller grains in low metallicity systems 
(ISandstrom et al. 2012h. since such dust g rains are more easily 
destroyed (e.g.. IWaxman & DraindbOOO ). However, this is an 
unlikely scenario for QSO DLAs, where sources of dust destruc- 
tion (shocks or radiation due to SNe or GRBs) are limited. Thus, 
the lower fraction of dust in the low-metallicity DLAs could be 
due to a lower efficiency of dust-production via grain growth and 
AGB stars, with little SN contribution. 



5.4. Dust extinction from the host galaxy ISM 

Given that the measurement of the dust extinction Ay is inte- 
grated along the whole line of sight, it is not trivial to assess 
whether the intrinsic Ay of the host galaxy is mostly associated 
with the bulk of the ISM. In principle there could be a contri- 
bution from intervening systems (although they s how typically 
weaker absorption than GRB-DLAs themselves; IVergani et al.l 
I2OO9I) or from a dense and dusty star-forming region associated 
with the GRB. On the other hand, we used the metal absorption 
lines to infer the dust column in the host galaxy ISM, specifically 
the column density of Fe in dust A^(Fe)dust (Eq.|9]l. It is insightful 
to compare Ay and A^(Fe)dust because they are two entirely in- 
dependent estimates of the dust content: Ay is derived from the 
continuum shape of the optical/UV afterglow, while A^(Fe)dust 
is derived from the optical absorption lines produced by met- 
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-1.0 -0.5 

log (Av/mag) 
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Fig. 6. The optical extinction Ay distribution with the column 
density of Fe in dust. The dashed line marks the column den- 
sity above which a signi ficant presence of molec ules has been 
observed in QSO-DLAs (iNoterdaeme et al.ll2008h . The symbols 
are the same as in Fig.|4] 



als at the redshift of the host galaxy ISM. Figure |6] shows that 
A^(Fe)dust steeply increases with Ay in GRB absorbers, at least up 
to Ay ~ 0.5. The significance of this correlation is confirmed by 
the Spearman's rank correlation test, as shown in Table |2] This 
strongly supports the general idea that Ay is mostly produced in 
the bulk o f the host galaxy ISM . A similar conclusion w as also 
reached bv lKruhler etaP fcoilh and lZafar et aH (1201 lab . 

GRB 070802 stands out from the observed trend with the 
highest Ay measurement. We note that this is one of the few 
GRBs where a Galactic -like 2175 A bu mp was observed im- 
printed on the opt ical afterglow spectrum (lElfasdottir et al.l2009t 
IZafar et all 1201 2h . The high extinction, well above the general 
trend of Fig. |6] might strengthen the hypothesis that in this case 
the extinction has been produced in two dif ferent regions with 
two different types of grains, as suggested bv lZafar et al.l(l2012l) . 
one of them not representative of the bulk of the host galaxy. 
One alternative explanation is that the distribution of A^(Fe)dust 
ffattens out, as if some saturation process prevented iron to be 
further accreted on dust grains above a certain threshold. 

In general, one obvious candidate for an additional region 
other than the host ISM to contribute to the Ay is a dense and 
dusty environment close to the GRB, where no low ionization 
lines can be produced. However, dust-destruction fronts typi- 
cally move faster and reach greater distances than the ionization 
fronts produced by GR Bs, making the dust destruction more ef- 
ffcient than ionization dDraine & Haoll2002l) . Thus, it seems un- 
likely to observe regions with dust but no low-ionization gas. 



5.5. Molecular content 

Give n the GRB association with massive s tar-forming regions 
(e.g. IWooslevllT99l iHiorth & Bloomll2012l) . we expect to ob- 
serve the signature of molecular clouds in some optical after- 
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glow spectra. In particular, dusty systems may bear the condi- 
tions for molecular formation, since dust grains can be efficient 
catalyst to aid the formation of mo lecules on their surfaces, such 
as H2 (e.g.. |Pirronello et"anil997l) . However, molecular absorp- 
tion bands are typically weak and notoriously difficult to ob- 
serve in optical/near-UV spectra of extragalactic sources (e.g., 
iNoterdaeme et al.ll2008h . In GRB afterglows CO was detected 
only once to da te in the low-resolutio n Keck/LRIS spectrum 
of GRB 080607 ( iProchaska et al.ll2009h . while H2 has been re- 
cently observed in a couple of mid-resolution VLT/X-shooter 
afterglow spectra (Kriihler et al., in preparation; Friis et al., in 
preparation). The dashed line in Fig. |6] shows the column den- 
sity log A^(Fe)dust - 14.7 above which molecules have be en ob- 
served in ~ 40% of QSO-DLAs ( INoterdaeme et al.ll2008h . Most 
GRB-DLAs lie above this column and thus are promising sites to 
detect molecules. 

While the presence of molecules is expected in GRB envi- 
ronments, the hard radiation of the bursts can be effective in de- 
stroyi ng dust and molecules in molecular clouds (iDraine &Had 
I2OO2I) . However, the currently low detection rate of molecules 
observed in GRBs may be an observational effect, due to the 
typically low spectral resolution or low S/N of the majority of 
GRB optical/near-UV spectra. On the other hand, current high- 
resolution samples might be biased against high metallicity lines 
of sight (potentially rich in dust and molecules) because only 
bright a fterglows can curren tly be observed at high spectral res- 
olution dLedoux et alj|2009h . A large sample of mid-resolution 
(e.g., VLT/X-shooter) spectroscopic observation GRB-DLAs 
can potentially include a significant number of those systems 
with the most favourable conditions for molecular formation 
(i.e., high H i and relatively high metallicity, iNoterdae me et al.l 
I2OO8I) . with the sufficient S/N and spectral resolution to constrain 
the actual molecular content and the extent of dust destruction. 

6. Conclusions 

We derive the dust-to-metal ratios (DTTVl) in long-duration 
GRB- and QSO-DLAs, based on the dust-depletion pattern ob- 
served from low-ionization metal absorption lines in the opti- 
cal/UV afterglow spectra. In this way the dust-to-metal ratio is 
purely associated with the the bulk of the neutral ISM in the 
GRB- or QSO-DLA, which is not influenced by the GRB, rather 
than integrated along the total line of sight. As a comparison, we 
also derive the UTM of GRB absorbers using a different method 
based on dust-depletion as well, described in ISavagliol (1200 Ih . 
with fully consistent results. 

The dust-to-metal ratios in GRB-DLAs span -10-110% of 
the Galactic value and show similar properties to the QSO- 
DLA distribution, with GRBs extending the study of dust-to- 
metal ratios out to higher column and volume densities. Above 
Z ~ 0. 1 Zq metallicities, these results are partially consistent 
with the finding of a univ ersally constant Galacti c-like dust-to- 
metal ratio based on Ay dZafarfe Watsonll2013h . However, at 
lower metallicities, down to Z ~ 0.002 Zq, we find compelling 
evidence for low values of DTyVt for both GRB- and QSO- 
DLAs. We argue that our UTM estimates for low-metallicity 
DLAs are more reliable than those based on the dust extinction 
along the line of sight because low Ay values are more sensitive 
to the contribution of foreground systems. 

The UTTvi distribution increases with metallicity and metal 
column dens ity. A simi lar trend was recently observed by 
iBrinchmann et al.l (1201 3l) for super-solar metallicity SDSS 
galaxies. We observe a flattening of this trend above log N{Zn) ~ 
13, suggesting that the production/destruction of dust with re- 



spect to the metals stabilizes for relatively more evolved galax- 
ies. The low dust fraction of low-metallicity systems may be the 
result of inefficient dust production, e.g. for grain growth and 
AGB dust formation in the low metallicity regime (where the 
time-scale for grain growth becomes too long), with little contri- 
bution by SNe. 

A possibility that needs to be further investigated is the pres- 
ence of a population of low-£)77Vl GRB-DLAs having a smaller 
amount of dust for their metal content, possibly due to destruc- 
tion of dust grains in actively star-forming regions. 

At z < 1.7 all the DTM estimates for the GRB-DLAs in 
our sample show high (Galactic-like) values. If confirmed with 
a larger sample, this may indicate that at lower redshift most 
GRBs occur in galaxies with Z > O.IZ© that are richer in dust 
with respect to the overall host-galaxies population. This could 
possibly be due to the overall increase of metallicity with cos- 
mic evolution. Interestingly, the presence of some of these high 
DTM galaxies is observed out to z ~ 4, indicating that for these 
systems a fairly rapid (< 1 Gyr) production mechanism is re- 
quired, i.e. in cooled SN ejecta, but also in AGB winds and grain 
growth in a Z ~ Zq metallicity regime. We find no trends of the 
DTM with the host-galaxy mass, star-formation rate and spe- 
cific star-formation rate. 

The column density of Fe in dust A^(Fe)du.st - derived from the 
relative abundances observed in the host-galaxy ISM - steeply 
increases with the dust extinction Ay observed towards GRBs, 
at least up to Ay ~ 0.5. This confirms that Ay is mostly pro- 
duced by the bulk of the ISM neutral gas. At high Ay we ob- 
serve one outlier of this trend, GRB 070802, where either an 
additional contribution for the extinction must be included or 
the A^(Fe)du,st flattens out at high Ay- The majority of GRB- 
DLAs shows high values of A^(Fe)du,st (above log A^(Fe)dust 14.7, 
the minimum value for QSO-DLA with molecular detection 
INoterdaeme et al.l [2008h . making them promising sites for the 
detection of molecules. 
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Ijenkinsi (l2009h . The slopes of the Galactic depletion se quences 
are c/Fec - Ape = -1.285 and dZnc - Azn - -0.610 (Ijenkinsi 
2009, where the depletion formalism is expressed in terms of A, 
B and z coefficients). The offsets of the Galactic sequences (i.e., 
depletion at Ft = 0) are offsetpe = Bpe - Ape zpe - -0.95 and 
offsetzn ~ 0, so then 



Szn.G = dZnc ■ Ft 

^Fe G = dFec ■ ^'h. + offsetpe. 



(A. la) 
(A. lb) 



Given that (5pe,G = <5zn,G + [Fe/Zn]G (see Fig. 1, regardless of the 
offsets), then Eq.|A]can be rewritten as 

^FcG - [Fe/Zn]G = dZnc ■ F.. (A.2) 

Deriving F» from Eq. lA.lbl and substituting in Eq. provides 

c/Znc 

^Fe,G - [Fe/Zn]G = — (5pe,G " offsCtpe). 

From dm - 1 - 10*' ' (Eq.O it finally follows that 

dtm(G) = 1 — j^QClFf^/Znlc i'Feo-offsetFc fl'ZnG)/(<iFeG-rfZnG) _ Q g(J 

(A.3) 



Appendix A: Derivation of the Galactic 
dust-to-metal ratio 

We calculate dtm{G) assuming [Fe/Zn]c - -0.95, which 
is the average value of Galactic disk-t-halo lines of sight 
d Savage & SembachI Il996h with the -0.03 dex correction of 
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[Zn/H] = -1.63 ±0.16 
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Fig. A.l. The depletion pattern [X/H] observed in GRB absorbers. The solid lines represent the best-fit depletion patterns corre- 
sponding to a best-fit dust-to-metal ratio, environment type and metallicity, labelled on each panel and reported in Table [1] [Zn/H] 
are used as metallicity input for the fit, while log Z/Z© are the fitted metallicities in case the H i column density but not A^(Zn ii) is 
constrained. The dashed lines refer to those cases where a single environment type could not be constrained. All upper and lower 
limits, as well as any constraint on sulphur are not included in the fit, but only displayed for completeness. 
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Fig. A.2. Same as Fig. lA.ll but for those cases where A^(H i) is not constrained, so the relative abundances with respect to Zn, [X/Zn] 
are fitted instead of the absolute abundances. 
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Fig. A.3. The DTTA distribution with the Hi column density. The symbols are the same as in Fig.|3] 
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